Thin film growth by using random shape cluster deposition 
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The growth of a rough and porous thin surface by deposition of randomly shaped clusters with 
different sizes over an initially flat linear substrate is simulated, using Monte Carlo technique. 
Unlike the ordinary Random Deposition, our approach results in aggregation of clusters which 
produces a porous bulk with correlation along the surface and the surface saturation occurs in long 
enough deposition times. The scaling exponents; the growth, roughness, and dynamic exponents 
are calculated based on the time scale. Moreover, the porosity and its dependency to the time and 
clusters size are also calculated. We also study the influence of clusters size on the scaling exponent, 
as well as on the global porosity. 



I. INTRODUCTION 

Nowadays, the fabrication of the solid-state electronic 
devices is vastly dependent on the deposition technol- 
ogy which is used in their construction. The formation, 
growth and dynamics of surfaces can well be regarded as 
the main factor in device fabrication which has attracted 
strong attention, specially in past few decades. Further- 
more, the deposition of particles via different procedures 
has immense importance not only in a broad range of 
physics applications, but also in chemistry, biology, and 
engineering. Until today, several analytical and numer- 
ical approaches and approximation methods have been 
developed to investigate the generic properties of various 
surface growth processes. Since the nature of resulting 
surfaces are random and far-from-equilibrium, they are 
called non-equilibrium and rough surfaces 1 ' 2 . 

Regarding it as the simplest model of surface growth, 
Random Deposition (RD) 1 has a continuum growth 
equation which results in exact determination of scal- 
ing exponents. In this model, each particle is randomly 
dropped over a site of an initially flat reference surface 
and deposited on the top of the selected column. The 
growth of columns are independent of each other, and 
the surface is uncorrelated. Ballistic deposition (BD) 3 
is another method of surface growth which is a colloidal 
aggregation. In this model, a particle is released from a 
randomly chosen position above the surface, located at 
the distance larger than the maximum height of interface. 
Since the different sites of the surface are dependent on 
the heights of neighboring sites, a correlation develops 
along the surface. Hence, the height of the new particle 
will be equal to or larger than that of its neighbors. The 
correlation length is the typical distance over which the 
heights "know about" each other which grows with time 
until it reaches the size of the system when the surface 
becomes saturated 1 . There are also other surface growth 
models such as Random Deposition with Surface Relax- 
ation (RDSR) 4 and nonlinear discrete models such as 
the Eden growth (ED) model 5 , restricted solid-on-solid 
(RSOS) model 6 and body-centered solid-on-solid (BC- 
SOS) model 7 have been proposed to study the kinetic 
roughening of the growth processes which belong to dif- 



ferent universality classes. 

In order the simple surface growth equations to present 
more complicated physical systems, theoretical efforts 
have been carried out and different types of noise were 
added to the models to adapt the possible variety of 
forces which might complicate the dynamics of growth 
in each physical system. The continuum linear growth 
model defined by the Edward- Wilkinson (EW) equation 
is counted with the same universality class as that of 
the discrete models, such as RDSR model 8 . The non- 
linear Kardar-Parisi-Zhang (KPZ) equation, which is a 
generalization of the EW equation, explains the univer- 
sality class which includes the discrete growth models, 
such as the BD, RSOS, and BCSOS models 9 . In some 
approaches, the two or more different deposition models 
are combined to investigate the time evolution of rough- 
ness in real systems 10-12 . The other works have proposed 
a surface growth model by the random deposition of lin- 
ear particles with different or identical sizes in (1 + 1) and 
(2 + 1) dimensions. These particles are aggregated on a 
fiat surface with special deposition rules and a porous 
bulk is generated 13,14 . 

Due to the daily development of nano-electronic tech- 
nology, an increasing demand is appeared to fabricate 
the nano-dimension devices. Deposition of clusters in- 
stead of single particles due to their large size, are a suit- 
able choice for building the desirable structures. Further- 
more, as the clusters occupy more than one unit site of 
the surface it is possible to generate porous bulks which 
has plenty of applications in magnetic storage 15 , solar 
cells 16 , carbon nano-tubes and other means 17 ' 18 . In the 
recent years, several methods have also been proposed to 
simulate the growth of thin films by deposition of clusters 
among which the Molecular Dynamics (MD) simulations 
and Monte Carlo (MC) technique have been frequently 
used. Unlike the MD model which is only applicable to 
small systems, the MC outstandingly can be applied to 
systems with large numbers of atoms 19-22 . 

In the following work, we are presenting a general- 
ization of the particle deposition to enrich the physics 
of the simulation process and make it correspond more 
with the real deposition processes such as sputtering in 
which clusters of random shape and size are ejected from 
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a solid target material due to bombardment of the target 
by energetic particles and deposited over the substrate. 
Although we do not deal with bombardment in our ap- 
proach, we propose a method to produce such random 
shape and size clusters, ranging from a single particle 
to a bunch of attached particles, to make the simula- 
tion corresponds more with experimental reality, we are 
also interested to the study of the morphology of sur- 
faces which are formed by adding such randomly shaped 
clusters to an initially flat linear substrate, based on MC 
simulations and random deposition model. 

In contrast to the ordinary RD in which there is no 
correlation along the surface and the interface width in- 
creases infinitely without occurrence of saturation, our 
approach results in aggregation of clusters which pro- 
duces a porous bulk with correlation along the sur- 
face. Moreover, the surface saturation is reached in long 
enough deposition times. Consequently, the scaling ex- 
ponents; the growth, roughness, and dynamic exponents, 
are countable based on the time scale. In addition, the 
porosity and its dependency to the time and clusters size 
is calculated. We have also studied the influence of clus- 
ters size on the scaling exponent, as well as on the global 
porosity. 

The outline of the paper is as follows: In section II, 
the model, deposition rules, and the details involved in 
the Monte Carlo simulations are described. The results 
of numerical simulations for scaling exponents, interface 
width, and porosity as a function of the time evolution 
and cluster's sizes are presented in section III. The paper 
end with conclusion In section IV. 



II. MODEL AND SIMULATION METHOD FOR 
CLUSTER DEPOSITION 

In the deposition process, the clusters with different 
shapes and size are dropped randomly on an initially 
flat linear substrate. As deposition starts, the surface 
grows gradually. The width of the surface which is called 
the surface roughness w, is a consequence of the non- 
equilibrium growth conditions under which the surface is 
formed 



w(L,t) = 
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where h(i,t) is the height of the growing surface at hori- 
zontal position i and time t and L is the size of the sub- 
strate. h(t) is the horizontally averaged surface height at 
time t 
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The comprehensive scaling behavior of the rough sur- 
face is studied with the Family- Vicsek scaling ansatz 23 . 



□ 







N=l N=2 



N=3 



N=5 



N=S 



FIG. 1. Pictorial examples of clusters with different sizes and 
random shapes 



The scaling relation implies that the interface width, in- 
creases as a power of time, w(t) ~ t& , which (3 is the 
roughness scaling exponent. The surface width saturates 
at a value which increases as a power low of the system 
size, w(t) ~ t a , where a is the growth scaling exponent. 
These two exponents are related to each other by another 
scaling exponent which is called the dynamic scaling ex- 
ponent; z 
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w(L,t)=L a f(^), 
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The function f{jr) is expected to have an asymptotic 
form such that 



w(L,t) 
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In addition to roughness width, w, the correlation length 
is another parameter which characterizes the rough sur- 
face by measuring the rate of changes of roughness along 
the surface. As the RD model has an uncorrelated na- 
ture, the correlation length is zero and the roughness ex- 
ponent is infinite. Therefor, the resulting surface is not 
self-afnne and the saturation will not be reach. 

In the following, at first, we will describe the proce- 
dure of building and depositing the clusters with differ- 
ent shape and size, using Monte Carlo simulation. Then, 
we show that during the surface growth, three different 
behaviors appear as a function of time. These behaviors 
at the initial and intermediate times, are characterized 
by growth exponents f3\ and /?2, respectively. At large 
times, the surface width grows slowly and the saturation 
will be reached. The morphology of rough surface is then 
characterized by the value of the exponent a. In addi- 
tion, the variation of global porosity with time and size 
of the bulk formed by depositing clusters are studied. 

In our approach for constructing rough structures, in- 
stead of growing the surface by atom-atom aggregation 
on the substrate, we propose a cluster growth procedure 
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FIG. 2. Illustration of the growth formalism and geometry 
conditions of the cluster deposition 
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in which a random number of single particles with a unit 
height are attached together in a way that different shape 
clusters are formed. Each cluster is a 3 x N linear box 
where N is the cluster length and ranges from one to an 
arbitrary random number. In other way, N is the maxi- 
mum permitted size for the deposition. The middle row 
of the cluster is fully filled with unit particles and the 
top and bottom rows are occupied by random number 
of particles, D, (0 < D < N), which their positions are 
quite random (Fig. 1). 

As the deposition is initiated, the clusters are dropped 
over random positions of a finite one dimensional sub- 
strate. The middle point of the cluster would land over 
the randomly selected position of the substrate. Regard- 
ing the size of a cluster, its middle point is considered as 
the half of its length. Thereafter, by checking the vacant 
sites around the landing point which depends on the ge- 
ometric constraint of the previous deposited layer (the 
previous MC step) the lowest possible place in which the 
cluster can be fitted is determined the incident cluster 
permanently sticks there. 

To clarify the growth rules and geometry conditions of 
the surface, three situations for which the clusters might 
be added to the surface are illustrated in Fig. 2. The 
vacant locations across the surface which has no other 
chance of being filled on the next deposition step, produce 
a porous structure in the bulk. The porosity is defined 
as the proportion of the unoccupied part of the surface 
volume to the total surface volume 



P = 



Vtoi 



(5) 



where V to t = V e + V s , which V s and V e are the volume of 
the occupied and empty sites, respectively. As the clus- 
ters are larger than unit site of surface, a porous struc- 
ture is produced such that its porosity is dependent on 
the time evolution and size of the clusters 24 ' 25 . In the 
next section, we focus on studying this dependency and 
scaling exponent of the surface by cluster deposition. 



FIG. 3. Log-log plot of the roughness as a function of time 
for deposition with different cluster shapes. From bottom to 
top the lattice sizes are L= 128, 256, 512, 1024, 2048. 



III. DISCUSSION AND RESULTS 

We performed a MC simulations on a linear flat sub- 
strate with it's size ranging from L = 2 7 to 2 11 . The 
deposition process starts at time t = and the log-log 
curve of the time dependency of roughness for different 
values of substrate and clusters size is evaluated. Fur- 
thermore, the behavior of the surface width and scaling 
exponents are studied for two cases of different and iden- 
tical cluster's size. The time evolution of the porosity and 
its dependency to the cluster's size are also discussed in 
section III C. 



A. Clusters with different size 

Considering the clusters with different size, we in- 
vestigate the time evolution of the surface roughness 
and estimate the scaling exponents and determine the 
universality class of the model. The MC simulations 
were performed for the substrates with different sizes as 
L = 128,256,512, 1024 and 2048. The numerical results 
show that as time goes by, in spite of the RD model, 
the surface roughness in our model presents three differ- 
ent behaves. This happens due to the correlation among 
the columns which leads to the lateral growth of the sur- 
face. The roughness behaviors correspond to the initial, 
intermediate and long time scales. The unit of time is 
measured in Monte Carlo steps (MCs). 

Fig. 3 demonstrate the log-log plot of the surface width 
as a function of time for the deposition of different size 
clusters on the substrate. The results were obtained for 
the deposition of a mixture of clusters with different sizes, 
3 x TV, where N ranges from 1 to 8. In order to get re- 
liable results, the average sampling were performed over 
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FIG. 4. Porosity as a function of time for depositing clusters 
with different shape and size on a linear substrate of size 
L=1024 



FIG. 5. Porosity versus cluster size for a linear substrate of 
size L=1024 



10 3 different samples for L — 128 — 1024 and 500 for 
L = 2048. As it is illustrated in the figure, the behaviour 
of surface width at initial times which corresponds to 
is close to J for all substrate sizes. This value is 
close to that of the random deposition of particles with 
arbitrary correlation. At the intermediate times, the lat- 
eral correlations are appeared and the growth exponent 
02 shows a completely different behavior. During this 
time, the height-height correlation develops and 02 ex- 
hibits smaller value than 0i . At this times 02 is approx- 
imately 1/4 which in contrast to 0i, does not belong to 
the RD universality. Although it is close to the results for 
the RDSR model in (1+1) dimension which is described 
by EW equation, however it is possible that our model 
does not belong to any universality class. It is due to the 
fact that the exponents have to be in exact agreement 
with that of EW so that we could decisively conclude 
that they are in the same universality class. Eventu- 
ally at long enough times, the a exponent characterizes 
the roughness of the saturated surface. The roughness 
width w has a power law dependency to the a expo- 
nents with respect to L (i.e.,u> ~ L a ). In other words, 
surface grows in early time regimes and saturates to a 
finite value, which depends on the size L. The exponent 
a is estimated close to 1/2. The scaling exponents were 
also measured for averaging over 1000 independent MC 
runs for each system size and the dynamic exponent z 
obtained equal to 2.00 from a and 02- 



B. Clusters with identical size 

Following our investigations, the MC simulations were 
performed for 7 different cases for each of which the clus- 
ters of identical size(iV) but with different shapes were 



TABLE I. The growth, roughness and dynamics exponents 
for different clusters sizes 



N 


Pi 


Pi 


a 


z 


2 


0.495±0.001 


0.293±0.001 


0.52±0.02 


1.77±0.07 


4 


0.532±0.004 


0.256±0.001 


0.50±0.03 


1.95±0.11 


6 


0.545±0.005 


0.260±0.002 


0.54±0.04 


2.07±0.15 


8 


0.553±0.007 


0.268±0.001 


0.55±0.03 


2.05±0.11 


10 


0.557±0.006 


0.267±0.002 


0.51±0.04 


1.91±0.15 


12 


0.562±0.007 


0.289±0.003 


0.58±0.05 


2.00±0.17 


14 


0.580±0.004 


0.296±0.004 


0.62±0.03 


2.09±0.10 



deposited over the substrate (Table I). 

The results show that the surface width does not ex- 
hibit a noticeable dependency to the size of clusters, and 
the time evolution of the roughness behaves in a similar 
manner to that of clusters of different sizes. In this case 
also there are two different regimes at the initial and in- 
termediate times which correspond to the exponents 0\ 
and 02, respectively. The exponent a is also determined 
at the long enough times when the saturation is reached 

Table I demonstrate the results of each MC simulation 
for different cases of deposition of identical size clusters. 
The scaling exponents 0\, /?2, a, and z for a substrate 
with L= 1024 are presented in the table. Although the 
clusters have identical sizes, the exponents behave the 
same as the previous approach with different size clusters. 
Once again we claim that although the exponents are 
close to that of EW universality class but as they are 
not in exact agreement, our model may not belong to 
any universality class. For example in the case N=4, 
the scaling exponents explain the EW universality class. 
However, it might be accidental. 
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FIG. 6. Image of the top most pieces of the surface which is 
obtained by our cluster deposition model. 



smaller than what is obtained from Eq. 6. Considering 
the N = \ case, after depositing a cluster on one side of 
the substrate, the probability of depositing another clus- 
ter exactly adjacent to it is very small (j-). Hence, the 
second cluster is mostly probable to land over the previ- 
ous one. As the height of each cluster in our model is 3, 
the total height of the surface at the end of the growth 
process is h = 3M. Inserting h in Eq. 6, the porosity 
reads P = |. For the case N = j , the second cluster sits 
in the vicinity of the first cluster with probability | and 
over it with | . The prosity is | for the former and | for 
the later and the total porosity is obtained as P = 0.71. 
Therefor, as the cluster size is decreased from -, the 
porosity is increased. The increase of porosity continues 
until the blending and mixturing of clusters revers the 
process and the porosity starts to decrease. 

Our investigations also shows that the porosity is in- 
dependent of the linear size of substrate. 



C. Porosity 

Deposition of particles of size 2 or larger would lead to 
the generation of porous structures 26-28 . Our model is 
also based on the deposition of different size and shape 
clusters and is expected to be porous. Fig. 6 illustrates 
a piece of the porous surface which is obtained from our 
deposition model. In the following, we investigate the 
porosity of our surface as a function of time and cluster 
size. 

Fig. 4 demonstrate the porosity of a surface of L = 
1024 with different size clusters as a function of time. Un- 
like the roughness width, porosity saturates at the early 
stages of deposition. The saturation occurs at 0.67. The 
dependency of global porosity of a surface of L = 1024 
to the cluster's size is depicted in Fig. 5. The results 
show that as the cluster size is increased, the porosity is 
increased either. And after reaching a maximum value at 
0.731, the porosity is decreased and tends to | in N = — . 
One can explain this behavior by considering the varia- 
tion of porosity in terms of the total height of the sur- 
face. If the total number of deposited clusters and the 
maximum surface height are denoted by M and h, re- 
spectively, the porosity can be approximated as: 



Due to consideration of the overhangs from the top- 
most layer, the maximum value for porosity is generally 



IV. CONCLUSION 

We have used the Monte Carlo simulations to study a 
surface growth model with deposition of the mixture of 
clusters with different sizes, 3 x N, where N is the clus- 
ter length and ranges randomly from one to a maximum 
value. 

Due to deposition of clusters with different sizes, a 
porous bulk is formed which reaches to a saturation 
regimes at the initial Monte Carlo steps (faster than the 
roughness width). The porosity depends on the cluster's 
size and increases as the size of clusters is increased until 
it reaches a maximum value, then decreases and tend to 
the value P = 0.67. 

Our results show that surface width presents three dif- 
ferent behavior as a function of time. At initial times, 
the behavior is close to an uncorrelated growth and the 
growth exponent is j3\ = 0.5. At intermediate times, the 
surface roughness grows more slowly with another growth 
exponent fii rangs 0.25 — 0.30. And eventually, at long 
times, surface reaches to a saturation regime which is 
characterized by the roughness exponent a rangs 0.5—0.6 
and as a result, the dynamical critical exponent rangs 
1.8 — 2.1. The results certify that, the present deposi- 
tion model does not belong to any universality class even 
though it's exponents are close to that of EW. 
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